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behaved, effective source for the corresponding waveform. 
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1. Introduction 

One of the main astrophysical targets of LISA pQ, a space-based interferometric 
gravitational wave detector, is the gravitational waves generated by the inspiral 
of compact objects into massive black holes. To extract physical information of 
such extreme mass ratio inspirals (EMRI), it is important to know the theoretical 
gravitational waveforms with sufficient accuracy. For the EMRI scenario, we use the 
black hole perturbation approach to compute waveforms. Here the compact object is 
approximated by a point particle orbiting a massive Kerr black hole. In order to obtain 
the precise theoretical gravitational waveforms we need to solve the self-force problem 
[21 EH H] and then for second order perturbations [51 H] . 

It has already been over two years that numerical relativity produced one of the 
most spectacular breakthroughs in science [HI [TJ [8], succeeding in solving the two body 
problem in general relativity after decades of effort. Among the notable set of results, 
the discovery [9] of large recoil velocities (up to 4000 km/s) [10] stands out. Notably in 
Ref. [9] a generic binary black hole case was treated, with unequal spins and unequal 
masses (mass ratio 1/2). It has been hard to deal with extreme mass ratios and recent 
computations limit to myjmi = 3/8 [TT] for spinning holes and up to nearly 1/4 for 
nonspinning holes [12J. It is foreseeable that soon mass ratios of nearly 1/10 can be 
achieved in this full numerical simulations making possible a comparison with the semi- 
analytic approach of the self-force. 

The self force problem was formally resolved over ten years ago [T3], [H], but 
its implementation in explicit computations proved lengthy and difficult. This have 
been reformulated in a more elaborated way by Detweiler and Whiting [15]. The first 
corrections to the trajectory of an EMRI was computed in [16j for a headon collision 
using the Regge- Wheeler gauge and the ^-function regularization. Those results were 
later confirmed using the standard formalism [T71 [181 IS] i n the Lorenz gauge [20J. 
In order to compute the generic orbit corrections around a Schwarzschild black hole, 
Barack and Lousto [21] have approached the problem of directly solving the linearized 
Einstein equations in the Lorenz gauge instead of the Regge- Wheeler and Zerilli wave 
equations [22]. To do that, one needs to be able to develop an accurate algorithm to 
integrate ten coupled wave-like equations with sources proportional to a Dirac's delta. 
This algorithm has been developed [231 121] for the (l+l)-wave equation resulting from 
the tensor harmonic decomposition of perturbations in the Schwarzschild background, 
but not for the (2+l)-equation resulting from perturbations of a Kerr (spinning) black 
hole. This is the subject of the current work. 

In this paper, we focus on one aspect of the self-force problem, specifically to derive 
the full, bare or retarded field of a point source. Therefore, we do not treat the local 
analysis of the field around the particle's location because the retarded field is global. 
Here, instead of studying the Teukolsky equation for the curvature perturbations ipi, as 
a first step, we consider the Klein-Gordon equation in the Schwarzschild spacetime, but 
do not decompose it into spherical harmonics, in order to model perturbations like in 
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the more generic Kerr background. Recently, introducing a thin worldtube surrounding 
the worldline of a point particle, Barack and Golbourn [25] have discussed this equation 
in (2+l)-dimensions as derived by the mode decomposition in the azimuthal direction. 
A different treatment is proposed here to deal with this problem globally. There is also 
a method to approximate a Dirac's 5-function by narrow Gaussian [26J. It is, however, 
difficult to ascertain the error introduced by smearing the particle and if this is accurate 
enough for self force computations. 

Once we obtain the retarded field, each azimuthal mode of the self-force on the 
particle can be calculated and is finite at the particle location. But the summation 
over all azimuthal modes diverges. Therefore, we need some regularization to derive 
the regularized self-force. At this stage, it is necessary to discuss the local analysis 
of the field or self-force in the derivation of the singular part. The regularized self- 
force includes two parts, i.e., the conservative part and the dissipative part [27] . To 
obtain the conservative part of the self-force, we need the regularization, while it is not 
necessary for the dissipative part which is derived by using a radiative Green's function 
|28j . Recently, Barack, Golbourn and Sago formulated a new scheme to construct the 
regularized self force directly from the azimuthal modes of the field in [29]. 

The paper is organized as follows. In section [21 we discuss the (2+l)-dimensional 
Klein-Gordon differential equation with a 2-dimensional ^-function as a source. To 
remove the 5-function, we introduce a new wave-function. This formulation is done in 
the case of general orbits in the Schwarzschild background. In section |3l we apply the 
formulation given in section [2] to the case of circular orbits. Here, we obtain a global 
effective source which is well behaved everywhere. To do so, we also treat boundary 
behaviors both near the black hole horizon and at spatial infinity. In section HI we 
summarize the results of this paper and discuss its applications. Some details of the 
calculations are given in the appendices. Throughout this paper, we use units in which 
c = G = l. 

2. Formulation 

When we calculate the (2+l)-dimensional equation derived from the 4-dimensional 
Klein-Gordon equation by the azimuthal mode decomposition, the resulting equation 
is not exactly same as the (2+l)-dimensional wave equation. In our formulation, it is 
important to derive a differential operator which is the (2+l)-dimensional d'Alambertian 
of the flat spacetime. By transforming the scalar field, we can obtain an equation which 
includes the flat (2+l)-dimensional d'Alambertian and a remainder as in (Q below. 
Then, we remove the 2-dimensional ^-function in the source term by using the Green's 
function method. 

In order to obtain the flat d'Alambertian, we consider the Schwarzschild metric in 
the isotropic coordinates, 

ds 2 = - f 2 P p l%y *? + (l + ^) 4 W + P 2 (d0 2 + sin 2 6dtf)] . (1) 
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This radial coordinate is related to that of the Schwarzschild, r, 



P 



r - M + (r 2 - 2Mr)^ 2 



(2) 



In the above coordinates, the Klein-Gordon equation with a point source reads 



(2, + M) 



16p 4 



(2p-M) 2 * ' (2/; • \lf 
16p 2 



d 2 + 



128/ 



(2p — M) (2p + M) 



+ 



(2p + M) z 



9^ + cot 6d e + 



sin 2 # 



x ifj{t,p,6,(f)) 



-q 



dr 



64p 4 5(t - ^(r))5(p - p 2 (r))5(0 - # 2 (t))<5(0 - 2 (r)) 



(3) 



(2p-M)(2p + M) 5 sin# 

where ■?/> and g denote a scalar field and a scalar charge, respectively. z a (r) is 
the particle's trajectory with a proper time r. Here, we use the azimuthal mode 
decomposition, 



ip(t,r,9,(j)) = ^ ^m(t,P,0)exp{ 



imc 



(4) 



and then obtain the wave equation for each m mode as 



(2p + M)' 



-d 2 
at 



16p 4 



{2p-Mf~ L ' (2p + M) 4 " 
16p 2 



128p 5 



(2p-M)(2p + M) 



5 d P 



+ 



(2 p + M) 
-Q 



di + cotOde- 



m 



sin 2 fl 



ipm(t,p,9) 



°° 64p 4 £(t - t z (r))£(p - Pz (r))5(6 - 
T (2p-M)(2p+M) 5 sinfl 
For the above equation, we transform the field ip m as 

, P 



exp[-2m0 2 (r)] . 



1/2 



1pm(t,P,9) 



Xm{t,P,0) 



(5) 



(6) 



_(2p + M) (2p-M) sinfl_ 
Above, we have set the radial and angular factors to obtain the spatial part of the flat 
(2+l)-dimensional d'Alambertian. Then, Xm satisfies the following equation 



1 (2p + M) { 
16(2p-M) 2 p 



2„4 * 



d 2 + d 2 + -d, 



+ - 



m 2 — 



l(4p 2 + M 2 ) 2 -16pMcos 2 fl 



Xm(t,p,6) 



sm 2 6 \" 4 (2p + M) 2 (2p-M) 2 
- <fr 2f( t - t ,(r))fO,-p,(r) W - fl,(r)) exp| _, m ^ (r)| 



[(2p - M)(2p + M)p sin 0] 1/2 

Note that from the relation z/> TO ~ Xm/p 1 ^ 2 , the source term for % m becomes p 1 / 2 times 
that for ^> m . This fact will be used in section [31 Next, we define a new time coordinate 



T 



1 dt A{2p z {t)-M)p z {tf 



{2p z {t)+Mf ' (8) 
where p 2 is obtained by solving the geodesic equation and we use the fact that the proper 
time r and the argument t in p 2 are also related by the geodesic equation. From this, 
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we derive an equation which can be divided into the (2+l)-dimensional d'Alambertian 
of the flat case IZ]( 2+1 ) and a remainder. 



C m Xm (T, p, 9) = (Q( 2+1 ) + CZ m ) Xm(T, p, 
= S m (T, p, 9) , 
where the differential operators are given by 

□ (2+1) = -% + d 2 p + -d p + - 2 dl , 

P P 



(9) 



prem 
*~'m 



(2p z (T)-M) 2 p z (T)\2p + My 



d 2 



(2p z (T) + M)6(2p-M)V 
2(4p,(T) - M){2p z {T) - M)(2p + Mfp z {TfM ( dp z (T) 



(2p,(T) + M)7(2p-M)V 
2 1 (4p 2 + M 2 ) 2 - 16pM cos 2 9 



dT 



p 2 sin 2 8 



m 



4 (2p + M) 2 (2p — M) 2 



(10) 



and the source term is shown to be 

2 5(t(T)-t z (T))S(p- Pz (r))5(e-e z (T)) 

dr -j- exp 



-im<p z (r)}- (11) 



'rem 



S m (T, p, 9) = —q 

\(2p- M)(2p + M)psm6Y 
Here, it is noted that there is no p and 6 derivatives in C r 
derivatives are included in □^ 2+1 ^ ) . 

To remove the 5-function in the source term, we set 

Xm (T, p, 9) = xl(T, P, 9) + X r r(T, p, 9) , 
where we define the new functions, Xm an d Xm™ as calculated from 
D^ X s m (T,p,9)=S m (T,p,9), 
C m xZ m (T, P, 9) = -CZ m X S m (T, P, 9) 

= S^\T,p,9). (14) 

The effective source contains no 5-function [j. Note that C^™ 1 includes a second- 

order derivative. But, since the factor of d\ is zero at the particle location, the singular 



of (JTOj) because all p and 9 

(12) 
(13) 



behavior of the effective source S, 



(e//) 



weakens. 



The derivation of the singular field Xm can ^ e performed through the Green's 
function of the (2+l)-dimensional flat case, 

1 



G(T,x;T',x') 



1 



1211/2 



9{{T~T') 



x'|), (15) 



2tt [(t - T') 2 - |x - , j 
where 9 denotes the Heaviside step function, x is a 2-dimensional spatial vector and the 
spatial distance is defined by 



p 2 + p' 2 — 2 p p cos( 



9') 



1/2 



(16) 



| This decomposition of Xm does not have any physical-meaning, i.e., Xm * s n °t identihed as the singular 
part to be removed in the self-force calculation. The physically and mathematically-meaningful singular 
part of the field which is called as the S'-part, have been discussed in [15]. Recently, Vega and Detweiler 
30| have discussed a new method to derive the retarded field by using a specific approximation to the 
S'-part. 
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Some detail on the above Green's function is discussion in section Appendix A Next, 
Xm is calculated by the following integral 



X S m (T,p,9) 



dT'p'dp'dO 1 G(T, x; T', x')5 m (T', p', 6') . 



(17) 



Here, it should be noted that the Green's function in the above multiple integration is 
given by analytically Therefore, it is easy to perform the integration in ( fT7l) even for 
general orbits. 

As for the remaining field Xm") we use numerical integrations. However, the 
effective source Sm^ is not amenable for direct numerical integrations due to its non- 
continuous behavior at the particle location. To explicitly obtain a source term well 
behaved everywhere, we will apply the formulation discussed above to the case of a 
particle on a circular orbit. 



3. Circular Orbit Case 



We consider a particle in a circular orbit given by 

z a (r) = {uV, r , |, u^r} , (18; 

where r denotes the orbital radius in Schwarzschild coordinates. The four velocity u c 
is written by 



u 



r - 3M 



1/2 



w 



M 



rg(r - 3M) 



1/2 



(19) 



The relationship between the new time coordinate T and the Schwarzschild time t 
can be obtained analytically 



where 



_ 4 (2p - M)pl 
(2 Po + M)3 ' 

r - M + (rg - 2Mr ) 1/2 
Po 7s 



(20) 



(21) 



Note that in general, for non circular orbits, we need a numerical integration to derive 
this relationship. 

In order to calculate the singular field, the Green's function in (JT5i) is rewritten in 
terms of t as 

G(T(t),x;T(t')X ' ' 



2tt 



16(2p,)-Af)2pg /, _ ,,x 2 _ I _ ,i 

(2p +M)6 y L b > \* A I 



1/2 



xfl |'*(fr-*)f8 (t _ 



x — X 



(22) 



(2p + Mf 

In the following, we will discuss the m = and m^fl modes separately. 

The effective source for a final regularized function Xm 9 must go like 0(p~ 2 ) for 
p — ► oo in the case of the m = mode and 0(p -3 / 2 ) for the m/0 mode. The reason 
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is the following. First, for p — > oo, there are only p n / 2 (n: integer) and In p terms in 
the effective source of this paper. And, in practice, we will calculate the regularized 
function of the original field ip m numerically. This means that the source term which 
we will use in the numerical calculation, becomes the factor ~ 1/p 1 ^ 2 times the source 
for Xm 9 ■ Thus, since we integrate the second order differential equation of (jHj), we need 
the above integrability conditions. In a similar manner, at the black hole horizon, i.e., 
p — > M/2, the source for i])™ 9 should vanish, i.e., the behavior of the source for Xm 9 
should be a power of (p — M/2) greater than 1/2 because the source for if) m is the factor 
~ (p — M/2) -1 ' 2 times the source for Xm 9 - 

3.1. The m = modes 

3.1.1. Singular field In the Schwarzschild background, the source term in ( ITTi) is time- 
independent, 

S (t, r, 9) = -q-— — 1 — -^5(p - p )5(6 - vr/2) . (23) 

M i [(2po-M)(2po + M)p ] 1/2 

Therefore, the solution x'o {t, t, 9) is obtained from the 2-dimensional Poisson equation 

as 



x S( tr 0) = _± h/£l — m(lx-xJ), (24) 

° l } 2 7 r nM(2 p + M)(2p -M)] 1 / 2 U Z ' J ' 1 1 

where the spatial difference is given by 

|x-x z | = (p 2 + Po- 2 ppo sin 6) 1/2 , (25) 

In practice, we must normalize the variables inside the logarithmic term in ( j24l) . We 
may take |x — x z | — > |x — x z |/M, but we will ignore the factor 1/M in the following. 

3.1.2. Local behavior The effective source for the m = mode, S^^\ is calculated by 
using the solution, Xoi as 



u* [(2p + M)(2p -M)] 1/2 P 2 sm 2 9 
(4p 2 + M 2 ) 2 -16p 2 M 2 cos 2 ^ _ 

(2p + M) 2 (2p-M) 2 ln(|X Xz|j - W 
This effective source is shown in figure [TJfor the innermost stable circular orbit (r = 6M) 
case, where q = 1, M = 1 and 6 = ty/2%. There is a singular behavior at the particle 
location. In order to perform the numerical integration with higher accuracy, it is 
convenient to regularize the source term to be at least C° at the particle location. 

In order to obtain the source S™ 9,1 which is regular at the particle location, we 
choose the regularization function Xo em,S as 

rem,s ( , m J_ pT (2 p - M) 3 ((4p 2 + M 2 ) 2 - 16p 2 M 2 COS 2 9) 
X ° {,P,) 16vr n*(2po + M) 5 / 2 (2po-M) n / 2 p9 

x |x-x z | 2 ln(|x-x z |) . (27) 
§ When we show figures for the m ^ mode, we also set t = and use the real part of sources. 
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Figure 1. Plot for the m = mode of S m with respect to p around the particle 
location. S { eff \ S r egJ , and So are shown by the solid red, dashed green, and 
dash-dotted black, respectively. The point particle is located at po ~ 4.95 (r = 6). 



This regularization function is chosen such that the additional source which arises from 
this function behaves well both for large p and at the horizon. Here, the factors of the 
power of 1/p and (2p — M) in the above equation give this good behavior. (We will 
discuss the m/fl case by using the same treatment for the regularization function.) 

Although this regularization function (and the other regularization functions 
discussed later) is clearly not unique, the sum of the regularization functions and the 
regularized function, i.e., the (original) retarded function has the physical-meaning and 
is unique. This is because the source for the regularized function changes by differences 
of the regularization functions. Therefore, one may construct any appropriate source 
by which we can derive a regularized function. 

Then, the regularized source at the particle location which is the source for 
Xo em ~ Xq 1 " 1 ' j is derived as 

S^^p^O) = S { eff \t,p,6) - C oX T' S ^P,0) • (28) 

We show the above source as the dashed green curve in figure HJ This regular source 
Sq £9,1 behaves as "xln \x\ for x — > 0", i.e., is C° around the particle location, and behaves 
as 0(l/p 2 ) for large p (See figure [2j). However, it diverges as 0(l/(p — M/2) 2 ) at the 
horizon (p = M/2) (See figure 0) . Therefore, we need one more regularization at the 
horizon. 



3.1.3. Behavior at the horizon Introducing the regularization function, 
X0l ' P ' ^ 2 7 r^ [(2po _ M )(2p + M)] 1 / 2 p 2 
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120 



Figure 2. Plot for the m = mode of S m with respect to p at large distance. S ' 



(e//) 



Sq 9 ' , and 5*0 9 ' J are shown by the solid red, dashed green, and dash-dotted black, 
respectively. 
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Figure 3. Plot for the m = mode of S m with respect to p near the black hole 
horizon. Sq*', S^ 9 ' 1 , and S^ 9 '* are shown by the solid red, dashed green, and 
dash-dotted black, respectively. The location of the horizon is p = 0.5. 



x - M 2 \nF(9) - -^—< ]rl L ' u ' 1 2 
4 V ; F{6) 



M 2 In F{6) - 4 p 2 Q In F(0) + 4 p sin M In F{6) 



-M 2 + 2 p sin M) - y ^ - (-M 4 In F(0) - 24 M 2 p 2 In F(0) 
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+8 M 3 p sin 9 In F{9) - 16 p 4 In F(9) + 32 p 3 sin 9 M In F(0) 
+ 16 M 2 p 2 cos 2 In F{9) — 3 M 4 — 44 M 2 p 2 , + 20 M 3 p sin 

+32 p 3 sin 9 M + 24 M 2 p 2 cos 2 s ) (p ~ y ) / ( m4 + 24 M 2 p 2 



-8M 3 p sin# + 16p 4 -32p 3 sin#M-16M 



2 Pg cos 2 



where 



F{9) = ^(M 2 + 4p 2 -4p sin#M), 

we derive the following source for the regularized function, Xo 
S r e9j (t, p, 9) = S r e9,I (t, p, 9) - C oX h (t, p, 9) 



rcy 



Xo 



rem,S 



Xo 



(29) 



(30) 

X h o- 
(31) 



Here, any bad behavior at infinity does not arise from the above regularization function 
because we have introduced the factor 1/p 2 in ( 1291) . This regularization function have 
been derived by using the Taylor expansion around the horizon, p = M/2 (This 
treatment is also used to derive Xm for m ^ in ( IB. II) .) We find that the regularized 
source is 0(p — M/2) as shown by the dash-dotted black curve in figure [Q, [2] and [3] and 
remains of 0(p -2 ) for large p. This completes the results for the effective source term 
of the m = mode to be used in numerical calculations. 



3.2. The m^O modes 

3.2.1. Singular field From f fTTl) . the singular field for the m^O mode is derived as 

(2p - M) 1 / 2 po /2 exp(-imu^T) 



X b m (t,r,6) 



2tt 



dr- 



(2p + My ft 



16(2pp-M)^ r t V; 

(2 P0 +A/)« ^ U T l 



X - X* 



1/2 



2tt 



' 4(2p -M)p 2 
X 1 (2p + M)3 {t UT) |X Xz 

K(2po-M)VVf 



x- 



«*(2p + M) 7 /2 
exp [—im(u^/ u l )T ] 



16(2p -Af)2p4 



x — x 7 



1/2 ' 



(32) 



(2p +Af) 6 

where the new variable for integration T , and the retarded time T ret are defined by 
To 



T -t ^ + M ) 3 lx x 



4(2p -M)p 2 
By introducing the following variable, 

_ 4 (2p - M)p 2 (t - Tp) 
(2p + M) 3 |x-x z | 



(33) 



(34) 
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2vrA M *[(2p + M)(2p -M)] 1/2 

exp{imn(2p + M) 3 |x - x z |T/[4 (2p - M)p 2 ]} 

X 



2 _ iW 2 



I 



(T 2 -l 

2 ^ exp(-zmOi) (1) / (2p + M) 



^3 



9 ttkH\ — -, — — -^-mfilx — x z | , (35) 

4 V[(2p + M)(2p -M)] 1/2 ° V4(2p -M)p2 1 z| 1 



where the angular frequency is given by 



n = H_, (36) 



and i^Q is the Hankel function of the first kind. The local behavior of the above solution 



near the particle location is 

X S m {t,r,6)~ x s m L (t, r ,6) 



2tt u t [( 2po + M)(2p - M)] 1/2 V4 (2p - M)p 



3.2.2. Local behavior When we write the singular field as 

(38) 

X„ is finite at the particle location. Then, the effective source in (114j) becomes 
S (eff) (t e) _ g 2 exp(-zmnt) / (2p + M) 3 

^ ( ' P ' )_ 2 7 r MM ( 2 , + M )(2p -M)] 1 / 2ln U(2po-M)pr fi|X ^ 



x I ^— I m 2 



1 (4p 2 + M 2 ) 2 - 16pM cos 2 9 



p 2 sm 2 e\ 4 (2p + M) 2 (2p-M) 2 
-(mfi) 2 f (2P0 + M)6 - ^ + M ) 6 ^ - £-v* (t r 0) (39) 

1 j V(2p -M) 2 p^ (2p-M)v/y m ml j ' 1 j 

Note that the third line of the above right hand side is at least C° at the location of the 
particle. Therefore, Sffi^ shown for the m = 1 mode by the solid red curve in figure HJ 
has a logarithmic divergence at the particle location. This behavior does not change for 
each m mode. 

To remove the logarithmic divergence in the source, we introduce 

x Po 19/2 (2p -M) 3 exp(^t) _ 32 

u* (2 po + M) 5/2 (2 po - M) 11/2 p 11 sin 2 # L 

+4 m 2 M 4 + 16 cos 2 6p 2 M 2 - (4p 2 + M 2 ) 2 ] . (40) 

Using this regularization function, we obtain a source S^ 9 ' 1 for the function Xm m ~Xm m ' S 

S^(t, p, 9) = S£M(t, p, 9) - C m X% m ' S (t, P, 0) ■ (41) 

The local behavior of S™ 9,1 , which is shown by the dashed green curve in figure HJ is of 
the form "xln |x| for x — > 0", i.e., C° at the particle location. 
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Figure 4. Plot for the m — 1 mode of S m with respect to p around the particle 
location. S[ eff \ S[ e3 ' J , S™ 3 * and S[ es ' f are shown by the solid red, dashed green, 
dotted yellow and dash-dotted black, respectively. The point particle is located at 
pa ~ 4.95 (ro = 6). The S^ 3 ' 1 , Sl eg ' h and S^ 9 '^ curves have an almost same behavior 
in this region. 



3.2.3. Boundary behavior We now focus on the behavior of the source term at the two 
boundaries, i.e., at the horizon of the large hole and spatial infinity. (See figure [5] and [6j) 
To regularize the source at the boundaries, we note that the source contribution from 
Xm" 1 ' 5 is we U behaved. This means that the ill behaviors of the source arise from x m - 
Therefore, it is convenient to use the asymptotic behavior of x m (and some correction 
factor) for regularization. 

For the regularization near the horizon, we use the regularization function x m given 
in section Appendix B Then, the source for the function Xm" 1 ~ Xm™' S — X m becomes 

S™ 9 > h (t, p, 9) = S^\t, p, 6) - £ mX h m (t, p, 9) . (42) 

This S^ 9,h is shown by the dotted yellow curve in figure [5] and vanishes as 0(p — M/2) 
at the horizon. But this source behaves as 0(p -1 / 2 ) for large p. (See figure [6j) To 
regularize it, we use the regularization function, 



1/2 



i q 



x p 



M 

T 



Pq 2 exp(-imQt) 

u t {2p + M) 2 (m ftp) 1/2 p 

i (2p + M) 3 mft (p 2 + pl 
4 



(P 2 + Po - 2 p p sin i 



exp 



2 po p sin i 



,1/2 



(2p -M)p2 
The final source for the regularized function Xm 9 

S™ 9j (t, p, 0) = S r ^ h (t, p, 0) - C mX ™(t } p, 9) . 



T 



(43) 



xZ m ~ xZ m ' S -X h m -xZ becomes 



(44) 
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Figure 5. Plot for the m — 1 mode of S m with respect to p near the black hole 
horizon. S{ e ^\ 5[ e9 ' 7 , S"[ e9 ' h and S[ eg ^ are shown by the solid red, dashed green, 
dotted yellow and dash-dotted black, respectively. The S l [ e9 ' 7 and S\ eg ' h curves have an 
almost same behavior except near the horizon. The location of the horizon is p = 0.5. 




Figure 6. Plot for the m — 1 mode of S rn with respect to p at large distance. S[ e , 
5[ e9,J , S\ e9,h and S^ 9 '-* are shown by the solid red, dashed green, dotted yellow and 
dash-dotted black, respectively. The S[ eff) , S{ e9j and S[ eg ' h curves have an almost 
same behavior in this region. 



This S^ 9 ^ is shown by the dash-dotted black curve in figure El and behaves like 
0(p~ 3 / 2 ) x (an oscillation factor with respect to p) for large p. Using this effective 
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source, we can calculate Xm 9 by numerical calculations. 
4. Discussion 

In this paper, we have formulated how to derive a global effective source, in replacement 
of a two-dimensional Dirac's delta, for the (2+l)-dimensional Klein-Gordon differential 
equation on a black hole background by using a transformation of the scalar field in ([6]) 
and a coordinate transformation with respect to time in (jSJ). Here, since we focus on 
the retarded field which is global, we do not use any local analysis and of the field, but 
have treated only the fields, i.e., x m i Xm™' 5 ) X m an d Xm defined globally. The above 
treatment of the regularization functions is the feature of this paper. On the other hand, 
Barack and Golbourn [25] have introduced a thin worldtube surrounding the worldline 
of a point particle, i.e., the local analysis to derive the retarded field. In this approach 
it is important to obtain results that are insensitive to the choice of the size of the 
world tube. On the other hand, our approach, determining a global effective source, is 
straightforward to use, once the regularization is done. 

As the application, in the case of circular orbit, we obtained the regularized effective 
source S^ 9 '? for the field Xm 9 - This source is C° at the location of the particle, and 
0(p — M/2) near the horizon. The behavior at infinity is 0(p~ 2 ) for the m = mode 
and 0(p -3 / 2 ) x (an oscillation factor with respect to p) for the m ^ modes, which 
allows straightforward numerical integration. 

In the case of general orbits, there is some additional difficulty. If it is possible to 
use the slow motion approximation and the eccentricity expansion for the bounded orbit 
cases, we can obtain x m analytically. Although we need to derive the singular field x m 
by numerical calculations in general, these include only the numerical integration for 
(jHJ) and (TTTT) . The regularization functions x m an d Xm f° r the field are extracted 
from the asymptotic behavior directly. This is the same method used in the circular 
case. About the regularization functions Xm™' S ) using ( TTTT) with the source which is 
evaluated from the asymptotic behavior of Sm , we can derive x™' s by a numerical 
integration because we need only the most singular part around the particle's location. 
Hence, there is no trouble to obtain the effective source for Xm 9 - 

When we consider the extension of this formulation to the Kerr background case, 
we can also extract a similar differential operator to that of (jHJ). In practice, we have 
the (2+l)-dimensional d'Alambertian of the flat spacetime in this case. The mode 
decomposition only in the azimuthal direction has already been used, therefore, the 
same treatment discussed in this paper is applicable. 

Finally, in the case of gravitational perturbations, we have ten field equations for 
the linear perturbation in the Lorenz gauge. (See [21].) We can also extract the (2+1)- 
dimensional d'Alambertian of the flat spacetime from them. The same treatment also 
holds to those equations. The method can also be used to deal with the Teukolsky 
differential equation [31] with the corresponding corrections for the presence of not only 
Dirac's delta but first and second derivatives of it as source terms. 
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Appendix A. About Green's function 

To obtain \m from ( p~3l) we use the Green's function method. Here we consider the 
Green's function in the Cartesian coordinates, i.e., x = rcosO and y = rsin# 

(_3 + dl + flj) G (t, x; t', x') = -S(t - t')6® ( X - x') , (A.l) 

where 5^ denotes the 2-dimensional ^-function. A solution of the above equation is 
usually calculated by using the Fourier transformation, 

1 f°° 

G(t, x; x') = — / duGufc t', x') exp(-iiut) . (A.2) 
2vr J-oo 

In the frequency domain, the Green's function satisfies the equation, 

(u 2 + d 2 x + d 2 y ) G w (x; t', x') = - exp(^t')£ (2) (x - x') . (A.3) 

The Green's function must become a homogeneous solution of the above equation for 
x 7^ x'. And, near the singularity, x — > x', the Green's function must behave as 
G w ~ — 1/ (2tt) In |x — x'|. Furthermore, we need to set a boundary condition. The out- 
going boundary condition is used here. From the above three conditions, the Green's 
function in the frequency domain is obtained as 

G„(x;C,x') = ^"Mx-x'DexpM'), (A.4) 



f i 

We go back to the time domain by using (1A.2j) . 



where Hq is the Hankel functions of the first kind. 



• POO 

G(t, x; t', x') = — / H^\u\x-x r \)exp[-iuj(t-t')]. (A.5) 
8vr ,/_ 00 

Note that H^\— \u\ |x— x'|) = — H^ 2 \\uj\ |x— x'|). Here, we use an integral representation 
of the Hankel functions 

rP)f~\ _ ~ 2i f°° exp(ixt) 



H^{x) = / dt 



7T j 1 {t 2 -iy/ 2 



r(2)/„\ 2 ^ /"°° Ji exp(-iit) 



7T 



H^(x) = - dt j-p > , (A.6) 



x {t 2 - 1)V2 



and then, the Green's function in the time domain is derived as 

G(i - x;t ' x '> = h [{t - t , Y X-A>r m " t0 " |x " x1) ' (A7) 

From the above Heaviside step function, this Green's function has support not only 
on the light cone, but also inside the light cone. This is mentioned as a failure of the 
Huygens principle [32] . This feature of the Green's function has been discussed in Refs. 
[33j [32l EH- The above Green's function is also derived from the direct integration of 
the (3+l)-dimensional Green's function with respect to the axial direction [33J. 
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The function Xm f° r ^ ne regularization near the horizon is so long that we summarize it 
in this appendix. 



v * (tnf)) = ± in v^ M4 

Xm[,P ' } 32 q ut[(2 Po -M)(2 Po + Mpz p 4 



exp(-imflt) 



where 
C (p, 6) 



x C (p, 9)H ( 



(i) 



-CM mi 



+C 2 (p, b)hP 



4(2p -M)p2 v ^ 
(2p + Mf 



(i) 



4 (2p - M)pl 



mny/F(6) 



(2p + M) 3 /— - 
4(2p -M)pr^ FW . 



(B.l) 



1 + 



8 

M 



M 



P - 



1 

64 



64 M 4 m 4 (2 p + M) 6 (M 4 + 12 M 2 p 2 cos 2 

-8 M 2 p 2 Q - 16 M cos 2 p 3 sin - 16 pi cos 2 + 16 p A )VL A 

+M 2 m 2 (29376 M 5 pjj cos 2 + 2448 M 7 pi cos 2 9 - 6274048 p^° 

+884 M s pi + 17 M 10 - 3264 M 6 p 4 cos 2 sin 9 - 43520 M 4 cos 2 9 p 6 sin 

-3128320 pi M 3 + 204 M 8 p 2 cos 2 - 52224 M p 9 cos 2 

-52224 M 2 p\ sin cos 2 9 + 1088 M 7 p 3 - 16320 M 5 pjj cos 2 9 sin 

+204 p M 9 - 17408 pj° cos 2 9 - 17408 M p 9 sin cos 2 

-65280 M 3 p 7 sin cos 2 9 - 272 M 7 p 3 cos 2 9 sin + 377984 p^ M 5 

-6343680 M 2 p\ cos 2 - 1540224 M 4 pg cos 2 + 6343680 M p 9 

+6287104 M 3 p 7 cos 2 9 + 11968 M 6 po cos 2 + 1629440 M 2 p 8 , 

+377984 M 4 p° - 102112 M 6 p 4 )fi 2 - 13824 p 4 (2p - M) 2 (M 4 



+ 16 M 2 p 2 cos 2 9 - 8 M' z p 2 + 16 p£) 



2 „2 



M 



[2p -Mf 



x (9 + 64 M m 2 fi 2 ) p 4 M 2 (M 4 + 16 M p 2 cos 2 — 8 M p Q + 16 p 4 ) 

d(p, 0) = i(2p + M) 3 fim(M 9 - 18p M 8 sin0 - 128 M 7 p 2 cos 2 + 144 p 2 M 7 

+448 pi M 6 sin cos 2 - 672 M 6 p 3 sin + 768 M 5 p 4 cos 4 

-2688 M 5 p 4 cos 2 + 2016 M 5 p 4 - 4032 pj> M 4 sin + 5376 p£ M 3 

+3584 pj> M 4 sin cos 2 — 512 pi M 4 sin cos 4 - 7168 p^ M 3 cos 2 

+2048 p£ M 3 cos 4 + 3072 p 7 M 2 sin cos 2 - 4608 p 7 M 2 sin 

AT 



+2304 Pq M - 2048 p* M cos 2 9 - 512 p 9 sin 0) ( p 
/ [(M 2 + 4 p 2 - 4 p sin M) 9/2 p 2 (-2 p + M) 

+ [64 M 2 m 2 (2 M 9 + 288 p 2 , M 7 + 7008 p^ M 4 sin cos 2 + 10752 p 6 M 3 
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- 14096 pi M 3 cos 2 9 - 8064 p 5 M 4 sin 9 + 4032 M 5 p 4 - 1024 p 9 sin 9 
-9216 p 7 M 2 sin 9 - 255 M 7 p 2 cos 2 9 + 4608 p 8 , M - 1344 M 6 p 3 sin 9 
+884 p 3 M 6 sin 9 cos 2 + 5952 p 7 M 2 sin cos 2 9 + 3904 p^ M 3 cos 4 
-4032 p 8 M cos 2 - 5316 M 5 p 4 cos 2 - 960 p 5 M 4 sin cos 4 9 
+1488 M 5 p 4 cos 4 9 - 36 p M 8 sin fl)ft 2 + 16 M 9 + 2368 p\ M 7 
+40960 Pq M + 33600 pi M 3 cos 4 5 + 59744 p 5 M 4 sin cos 2 9 
-80896 p 7 , M 2 sin 9 - 11200 M 6 p 3 sin - 121872 pg M 3 cos 2 
-44804 M 5 p 4 cos 2 + 52032 p 7 M 2 sin cos 2 9 - 35776 p 8 , M cos 2 5 
-9216 p 9 , sin 9 + 93184 pg M 3 - 8128 pg M 4 sin 9 cos 4 

„3 i\/r6 ■„ n „„„2 /) eonno „S ji,f4 •„ /i , i o/in« y\/f5 „4 „„„4 , 



+7348 pg M° sin 9 cos z 9 - 68992 pg M sin + 12496 M s p 4 cos 4 9 
-292 po M 8 sin 9 - 2095 M 7 p 2 cos 2 9 + 34048 M 5 po rn Q (2 p + M'f 
m \~ 

x I P 



11 " ' (-2po + M)Mp 2 (M 2 + 4p 2 -4p sin#M) 9/2 



2 , 

x (9 + 64 M 2 m 2 Q 2 ) 

C 2 (p, 9) = ^-(M 2 -4p sin^M + 4p 2 -4 cos 2 ^p^)fi 2 m 2 (2p + M) 6 (64M 2 m 2 fi 2 
64 

+!)(p-y) /[^)(-2p + M) 2 p 4 (9 + 64M 2 m 2 fi 2 )] (B.2) 
In the above equations, F(9) is the same as defined in fl30|) . 
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